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Abstract 

We calculate the spectrum of the scattered light from quantum degenerate 

atomic gases obeying Bose-Einstein statistics. The atoms are assumed to 

occupy two ground states which are optically coupled through a common 

excited state by two low intensity off-resonant light beams. In the presence 

of a Bose condensate in both ground states, the atoms may exhibit light 

induced oscillations between the two condensates analogous to the Josephson 

effect. The spectrum of the scattered light is calculated in the limit of a 

low oscillation frequency. In the spectrum we are able to observe qualitative 

features depending on the phase difference between the macroscopic wave 

functions of the two condensates. Thus, our optical scheme could possibly 

be used as an experimental realization of the spontaneous breakdown of the 

U(l) gauge symmetry in the Bose-Einstein condensation. 
42.50.Vk,03.75.Fi,05.30.Jp 
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The Bose-Einstein condensation (BEC) has finally been observed in a weakly interacting 
system with well-understood interactions by cooling and trapping alkali metal vapors |1] [| . 
In these first experiments on BEC the condensates were probed destructively with a short 
flash of laser light. Recently, Andrews et al. f|] have reported a non-destructive optical 
detection of a Bose condensate. This detection technique could possibly also be used to 
measure the spectrum of the scattered light from the condensate, when the driving light is 
detuned far from the resonance of the optical transitions f5^0. 

In the BEC phase transition the Bose gas should acquire nontrivial phase properties. In 
the traditional reasoning the condensate is given a macroscopic wave function which acts as a 
complex order parameter with an arbitrary but fixed phase ||. The selection of the phase in 
each experiment implicates the spontaneous breaking of the U(l) gauge symmetry. However, 
even though the condensate is taken to be in a number state with no phase whatsoever, the 
condensate behaves as if it had a phase fl9[-|i6|j . The phases of the condensates are expected 



to exhibit collapses and revivals due to collisions ||n|,|T 

Most of the theoretical studies on optical response of degenerate atomic gases p|-fr|,|TE|-|2"l 



are insensitive to the phase properties of the condensates. Javanainen PS[ and Imamoglu 



and Kennedy |2?J have proposed optical schemes to detect spontaneous breakdown of the 
gauge symmetry in BEC. In Ref. f25[ two condensates in different Zeeman states are confined 



in the same trap, and two phase coherent laser beams are used to drive Raman transitions 
between the condensates. Amplification of one of the light beams, because of the stimulated 



Raman scattering, would indicate phase dependent properties. In Ref. |26| nonlocal light 
scattering between two independent, spatially separated condensates were studied. 

In this paper we consider the spectrum of the scattered light from a quantum degenerate 
Bose-Einstein gas occupying two ground states with a common excited state. Following 
Ref. [^5j all the atoms are confined in the same trap. The possible Bose condensates are in 
two different Zeeman sublevels |6) = \g,m) and |c) = \g,m — 2). The state c is optically 
coupled to the electronically excited state \a) = \e,m — 1) by the driving field E^ 2 having 
a polarization a + and a dominant frequency f^. Similarly, the state b is coupled to a by 



the driving field E^i with a polarization cr_ and a dominant frequency We assume that 
the light fields E dl and E d2 propagate in the positive ^-direction and are detuned far from 
the resonances of the corresponding atomic transitions. The light fields are also assumed to 
be in the coherent states. In accordance with our previous work 0,0] we can immediately 
write down the Hamiltonian density 

H(r) = ^(r)# c . m .^(r) + ^ c {r){H c . m . + fwJ cb )^ c {v) + ^ a (r)(H c . m . + huj ab )^ a {r) + H F (r) 
^d 6a .E 1 (r)^ t (r)^ a (r)+h.c.) - (d ca • E 2 (r) ^(r)i(r) + h.c.) , (la) 



Hp — j H f (t) = hJ2 u q a\a q . (lb) 

The first three terms reflect the energies, internal and center-of-mass (cm.), of the atoms 
in the absence of electromagnetic fields. The frequencies for the optical transitions a — > b 
and a — > c are u ab and u ac (u c b = tOab — uj ac ), respectively. The last two terms are for the 
atom-light dipole interaction. The dipole matrix element for the atomic transition a — * b is 
given by d& a . The Hamiltonian density for the free electromagnetic field is Tip- If the light 
field is fully treated as a quantum mechanical field with its own dynamics, the Hamiltonian 
may also contain the polarization self-energy term and the dynamical light field should be 
the electric displacement instead of the electric field [^4],^7]]. However, if the detunings 
(Ai = Qi — Ld a b ~ A 2 = fi 2 — ^>ac) are sufficiently large [p5| , multiple scattering is negligible 



and, at first order in 1/Ai, these effects may be ignored. 

We define slowly varying field operators in the Heisenberg picture by ip a = e tnit tp a , 
E+ = e init Bf, E^ = e in2t E£, and ip c = e i{ni ~ Q2)t ip c . In the limit of large detuning, the 
excited state field operator ip a may be eliminated adiabatically: 

Mr) = ^- (d ab ■ E+(r) ^(r) + d QC • E+(r) &(r)) . (2) 

We insert Eq. (0) into the Hamiltonian (|2|) and keep only the terms of first order in 1/Ai. If 
the detuning from the two-photon resonance is given by S c b = Qi — fi 2 — ^cb, the Hamiltonian 
density now reads 



H = i)\H c , m ^ h + ft c {H c . m . - h5 cb )4> c + ft F - _ {d ba ■ Er dafe • E+ ^ b 

+ d ac • E+ d ca ■ E 2 - ftfc + (d ba ■ d ac • E+44 + h.c.) } . (3) 

Then, the electric fields may be solved and, according to Ref. ||, the scattered fields are 
given by 

KiW = J^-J d 3 r'K(d ba ; r - r') {d afe • E+(r) ^(r')^(r') + d ac • E+(r) ^(r')4(r')} , 

(4a) 

E+ 2 (r) = ^1- J rf 3 r'K(d ca ; r - r') {d ac ■ E+(r) ^(r')^c(r') + d ab • E+(r) ^(r')^(r')} . 

(4b) 

Here we have used the first Born approximation based on the assumption that the incoming 
fields dominate inside the sample as the multiple scattering is negligible. The kernel K(T>; r— 
r') is the familiar expression [^HJ of the positive-frequency component of the electric field 
from a monochromatic dipole with the complex amplitude T>, given that the dipole resides 
at r' and the field is observed at r. 

We assume the driving electric fields to be plane waves 

E+(r) = \£ x e^\ E+ (r) = \£ 2 e + e^ . (5) 

We insert these into the Hamiltonian (|3]) to obtain the Hamiltonian for the matter field 
dynamics 

Km = 4(H c .m. - n8 1 )ij h + &(H c . m . - H5 cb - hfo)#c + (h^l^e'^ + h.c.) , (6) 
where k 12 = n\ — k 2 is the wavevector difference of the incoming light fields. Following the 



notation in Ref. 



|£i| 2( ^ab x _ \£-2\ 2 d 2 ac _ £*£2d ab d ac 

5i "4^a7' 52 "4^a7' k ~ WK, ' (7) 

The dipole matrix element d ab contains the reduced dipole matrix element and the corre- 
sponding nonvanishing Clebsch-Gordan coefficient. To simplify the algebra, we assume k to 
be real. 



We assume a translationally invariant and non- interacting Bose gas as in Ref. ||25|| . 
The matter field operators are given by the familiar plane wave representations ^(rt) = 
y-i/2 ^ k e * k r bk(t) and ^(rt) = V~ 1 ^ 2 J2k e * k r Ck(t). In the absence of light, the cm. motion 
in both ground states satisfies the dispersion relation for the frequency e k = h\k\ 2 /2m. 

By assuming that K12 is very small we may neglect its contribution to the time evolution. 
The dynamics of the ground state annihilation operators Ck and 6k rnay then be easily solved 
from the Hamiltonian (||). With the definitions 28 = 5 c b — Si + 5 2 and Qr = (5 2 + k 2 ) 1 ^ 2 the 
solutions are given by 

= e^+s^t | 4(0) L 0S Q Rt + ^_ sin Q R A _ ^ &k _ Ki2 ( ) sinfi^j , (8a) 

b k (t) = e^-^ j& k (0) (cosSl R t - ^- sinfifltj - ^-5 k+Kl2 (0) sin^tj . (8b) 

Before the light is switched on, the atoms in the states b and c are assumed to be uncor- 
rected. The driving light fields induce a coupling between the two levels. In the presence 
of Bose condensates in the ground states, the coupling between the two condensates is 
analogous to the coherent tunneling of Cooper pairs in a Josephson junction P^ipS . 



According to Ref. 0, the spectrum is obtained by calculating correlation functions for the 
matter field operators. In particular, in the presence of the condensate the expectation values 
for single operators are nonvanishing (0 C (O)) = (N c /V)^ 2 e Wc and ("06 (0)) = (N b /V) 1 ^ 2 e lLPb . 
We have assumed that the condensates in the states c and b have the expectation values 
for the number operators N c and Nb, respectively. The condensate phases tp c and tp^ are 
random variables and are fixed in each experiment. 

From the equations (|j) and (|8]) we immediately see that the spectrum of the scattered 
light is time- dependent. However, with the measurement times short enough compared to the 
characteristic time scale of the changes in the spectrum, the spectrum may be considered 
steady. Evidently, we do not want too short measurement times which would wipe out 
qualitative features from the spectrum. According to Ref. 0, the most relevant choices 
for the characteristic frequency scales in the spectrum are the effective recoil frequencies 



= h(AKi) 2 /2m, where A/^ = fljfi/c — Kj is the change of the wave vector of the light 
field Ej upon scattering and n is a unit vector pointing into the direction of the scattered 
light. 

We present a set of assumptions based on the limit of a weak coupling between the 
ground states to obtain a stable spectrum during short measurement times. We assume 
that the oscillation frequency Q R is small, u Ri ^> Q R , or u Ri ^> 5 and oj Ri ^> k. The latter 
indicates AiUm 3> r R,\R,2- t where 'Kx = d a bSi/2h and IZ2 = d ac £ 2 /2h are the Rabi flopping 
frequencies for the electronic transitions b — > a and c — > a, respectively. The quantities u R i 
are expected to be in the neighborhood of atomic recoil frequencies, um ~ e R ~ 100 kHz as 
a rule of thumb ||. The detuning Ai may be chosen large keeping in mind that the intensity 
of the scattered light scales as 1/A 2 . If the detuning from the two photon resonance becomes 
small, the effective linewidth 7 of the transition c — > b may have an effect. However, it may 
be shown to be proportional to A{ 2 or smaller. Thus, we can also safely assume cum 7. 
The requirement, for the measurement times T to be short enough, sets the conditions 
kT < 1 or > TIZ1IZ2 and 8T <C 1. With these assumptions the time evolution for 
the atom operators from Eq. (^) with the short measurement times becomes very simple: 
c k (t) = 4(0) e i(5l - e k)* and b k (t) = 6 k (0) e i(<5l " e k)*. 

The two components of the scattered light E si Eq. fliaD , and E s2 Eq. (|4b|), have polar- 
izations cr_ and a + , respectively. By selecting the polarization we may detect these fields 
separately. Let us first consider the spectrum of the scattered field E sl . To detect only E sl , 
we need to select the polarization e = (\/2e_i —tan 9 e~ l ^e z ) / V tan 2 9 + 2 from the scattered 
light, where 9 is the scattering angle in the conventional spherical representation (the angle 
of the scattered light with respect to the positive z-axis). The polar angle <fi is the angle with 
respect to the x-axis in the xy-plane. This angle has an influence on scattering processes in 
which the phases of the two circularly polarized light beams are important. The calculation 
of the angular distribution of the scattered light with the given polarizations and atomic 
transitions is explained in Ref. |J. If we neglect the forward scattering, immediately after 
switching on the light the spectrum is given by 
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5 1 (r;«) = C7 1 (r)[5 1 (a;)+5i6(«)] l Ci(r) = 2fl "ffi 2A o • ( 9a ) 

16 (tan ft + zjTT^AieocrV^ 

«Si(u;) = 5i <J(a; + uj r1 ) N b (l + n bAK J +6 1 5(w- u m ) N b n' bAKi + 5 2 5(uj + u R2 ) N c (l + n' bAK2 ) 



+ 5 2 5{u - u R2 ) N b n' cAK2 + Ik cos(A^) 5{uj + u m ) 5(A Kl - Ak 2 ) yj N c N b {l + n' bAK1 ) , (9b) 

S lb (u)= 5i S(u - e k + e k _ AKl ) (1 + n' b>k _ AK1 ) 
k 

+ ^2 5(w - e k + e k _ AK2 ) < k (l + ^.k-AKa) • ( 9c ) 
k 

We have explicitly separated the contribution of the Bose condensates into S\(u). This 
term is nonvanishing only, if there is a condensate present at least in one of the states b or 
c. The term Si b (u) arises from scattering processes between the noncondensate atoms. The 
expectation values for the number operators for the states c and b in the cm. state k are 
n c k = (z _1 e /37i ^ k+Wcb ' ) — and nj, k = (z~ 1 e l3nek — l) -1 . The primes in these terms indicate 
that the condensate states are excluded to avoid double counting. For example, the notation 
in the last term in Eq. (|9b|) n' bAKl should be interpreted as Ak,\ ^ 0. The difference between 
the condensate phases is written as A(f — <p b ~ fc- We have used the same normalization in 
Eq. (g) as in Ref. §. 

The first term in Eq. ( |9"c| ) describes a scattering process in which an atom in the ground 
state b with the cm. state k scatters to the cm. state k — Aki still remaining in the state 
b. The delta function dictates the energy conservation, which coincides with the theory 



for Doppler velocimetry of atoms [§0J shifted by the effective recoil frequency ujri . The 
product of the occupation numbers indicates that the scattering is enhanced if the final state 
is already occupied. The first term in Eq. (^c|) and the first two terms in Eq. (|9"15|) coincides 
with the results given in Ref. 0, apart from the angular distribution of the scattered light 
and possibly the Clebsch-Gordan coefficients of the atomic transitions. In the second term 
in Eq. fl9"c|) an atom scatters from the ground state c to the ground state b, while the cm. 
state undergoes the change k — > k — An 2 - 

The last term in Eq. ( |9"b~D is purely a consequence of the broken gauge symmetry in 
BEC and it depends on the phase difference Atp between the two condensates in the states 



b and c. It varies from measurement to measurement, because Aip is essentially fixed in 
each experiment as a random number. The delta function for momenta indicates the exact 
conservation of the momentum of the two scattered photons for an ideal gas in a spatially 
homogeneous system. For a weakly interacting gas in a magnetic trap, the momenta of the 
scattered photons do not need to be exactly the same. 

The last term in Eq. (|9b|) is strongly peaked at uo = —uo R \. With An\ = An 2 and 
by setting N b = N c and 5\ = 5 2 = k the condensate part of the spectrum reduces to an 
especially simple form: 



= 25 1 5(u + u R ) N b (l + n' bAK )[l + cos (A<p)] + 5 1 5(lu - u R ) N b {n' bAn + n' cAK ) . (10) 



The effect of the phase difference between the condensates is clearly observed. The peak at 
u = —ujr reaches its maximum at A(f = and completely vanishes at Aip = tt. 

To detect only E S 2, we need to select from the scattered light the polarization e = 



the same as Si(ru) given by Eq. (|^) with the following obvious changes: c <-> b, Ak,i <-> Ak 2 , 
and Q2 — * 

If the detector is insensitive to polarization of the scattered light, the spectrum has 
contributions from <Si(r;o;), <S 2 (r;u;), and from the cross terms iS>i 2 (r;c<;) oc (E^E^) + 
(E^Ejg). These cross terms depend on the phases of the two circular polarizations and 
the scattered light does not have a rotational symmetry around the z-axis. Neglecting the 
forward scattering, the resulting spectrum, immediately after switching on the light, reads 



(y^e+i + tan#e^e z )/"\/tan 2 9 + 2. The scattering spectrum ^(rtu) obtained in this case is 



Stir; u) = C[(v) [S x {uj) + S lb (uj)} + C' 2 {v) [S 2 (cu) + S 2b (u)} + S 12 (tu) , 



(11a) 




n\dl h {2 - sin 2 9) 
647r 2 A 1 eoC 3 r 2 




njdl(2-sm 2 9) 
Q4ir 2 AieoC 3 r 2 



(lib) 



QlQld ab d ac sin 2 9 



cos 



(2(f) + Ay?) y/AyVc (5i 5(u - oj r1 ) n' bAni + 5 2 8(u - u R2 ) n' cA ^ 2 ) 



327r 2 Aie c 3 r 2 



(11c) 



8 



The spectra Si(u) and Sib(u) are given by Eq. (J9|) and the spectra £2(0;) and £25(0;) may 
be obtained from the same expressions with the previously explained changes. The new 



term (lie) also is purely a consequence of the phase difference between the two condensates. 
By setting Aki = A/t 2 , N b = N c , d ab = d ac , 5i = 62 = k, and the scattering angle 
9 = 7r/2 the total scattered spectrum may be expressed in the simple form: S(r; uo) = 
C[(r) [S(u) +S b (u)\, where 

S(u) = 26, 6(uj + u R ) N b (2 + n' bAK + n' cA J[l + cos (Ap)] 

+ 25 1 5(lu - u R ) N b (n' bAK + n' cA J [1 + cos (20 + Ay?)] , (12a) 

SftM = 2 51 J l - e k + e k-A K ) n' ak (1 + ^, k -A«) • ( 12b ) 

Here the summation over a and (3 denotes the summation over the ground states h and 
c. The phase difference between the condensates again has clearly observable effects on the 
spectrum. The two condensate peaks in Eq. ( |12a| ) oscillate as a function of the phase differ- 
ence between the two condensates. However, these oscillations depend on the observation 
point of the scattered light (the angle <fi in the xy-plane). The detection of the scattered 
light at several different points in space (with different values of 0) may provide a scheme 
to determine the phase difference between the condensates very accurately. 

The integration of Eq. (|12b|) is straightforward in the continuum limit and the results are 
given in Ref. 0. For simplicity, we choose uj c h = 0, so that the statistical distributions of the 
states b and c are identical. We give representative spectra of Eq. (|T2"|) in Fig. [I] with different 
values of Atp. We choose the detection point at 9 = it/ 2 and <fi = tc/2. To simplify further, 
we set ujr = ujd = [ksT '(Ak) 2 / m] 1 / 2 , where u>d is the effective Doppler width corresponding 
to the change of the wave vector Ak. For 87 Rb atoms this corresponds approximately to 
the temperature 270 nK. We remove the integrable singularities of Eqs. ( |1 2a| ) and ( |12b| ) by 
convoluting the spectra with a Gaussian whose variance is 0.01 ujr. It is assumed that a 
fraction 0.8 of the atoms in both ground states are in the Bose condensate. 

At Aif = the peak at u = — ujr has its maximum height. The contribution of the 
scattering processes into the condensates completely vanishes in the spectrum. The peak at 
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u) = ujr has its minimum height. At Aip = tt /2 the cosine terms in Eq. (|12a|) vanish. At 
Aip = 7i the contribution of the scattering processes out of the condensates vanishes in the 
spectrum. The peak at u = — cur reaches its minimum, while the peak at u = ojr is in its 
maximum. 

The spectra (^) and (|TTJ) were calculated immediately after switching the light on. The 
expectation value of the atom number in the two ground states oscillates with the frequency 
2Qr, see Eq. (||). The atom numbers of the condensates oscillate even if the number of 



atoms in each condensate is initially equal, analogous to the Josephson effect [P9| , P5|| . With 
5 ~ and — the expectation values for the number operators in the two condensates 
would approximately interchange in the time period of tt /2k. 

We have calculated the spectrum for an ideal, noninteracting gas. In a magnetic trap and 
in the case of a weakly interacting gas the momentum distribution is broadened, due to both 
the uncertainty principle and the interparticle interactions. However, the qualitative features 
in the spectrum are still well recognisable |||7]|. If we wanted to consider Bose condensates 
in arbitrary momentum states, we would evidently have to abandon the assumption of 
translational invariance in space and use some other eigenfunction basis than the plane 
wave basis for matter field operators. 

In conclusion, the spectrum of the scattered light from two Bose condensates in a trap 
may show an unambigious signature of the broken gauge symmetry in BEC. The phase differ- 
ence of the macroscopic wave functions for the condensates has clearly observable effects on 
the spectrum and its value may possibly be determined very accurately from the oscillating 
heights of the peaks. Recently, Myatt et al. have managed to produce two overlapping 
Bose condensates in two different angular momentum states of 87 Rb. The condensates were 
created using nearly lossless sympathetic cooling of one state via thermal contact with the 
other evaporatively cooled state. The two states \F — l,m = —1) and \F — 2,m — 2) 
cannot be coupled by two photon transitions, as was assumed in our calculations, but the 
generalization for the three photon case should be rather straightforward. 
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FIGURES 

FIG. 1. Spectra of light scattered from a Bose gas occupying two ground states with various 
values of the phase difference Aip between the two condensates. The scattered light is measured at 
6 = 7r/2 and (j> = tt/2. A fraction 0.8 of the atoms in both ground states are in the Bose condensate. 
For the solid line Aip = 0, for the dashed line Aip = tt/2, and for the dotted line Aip = n. The 
origin corresponds to the common laser frequency f^i = f^. The spectra are expressed in units 
of the effective recoil frequency lor, and the effective Doppler width is chosen as ur> = ojr. To 
regularize integrable divergences, all computed spectra are convolved with a Gaussian with the 
variance 0.01 ujr. 
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